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(57) ABSTRACT 

Method and system for monitoring and analyzing, in real 
time, variation with time of an aircraft flight parameter. A 
time-dependent recovery band, defined by first and second 
recovery band boundaries that are spaced apart at at least one 
time point, is constructed for a selected flight parameter and 
for a selected time recovery time interval length At(FP;rec). A 
flight parameter, having a value FP(t=t ) at a time t=t , is 
likely to be able to recover to a reference flight parameter 
value FP(t';ref), lying in a band of reference flight parameter 
values FP(t’;ref;CB), within a time interval given by 
t / ,^t'^t p At(FP;rec), if (or only if) the flight parameter value 
lies between the first and second recovery band boundary 
traces. 

20 Claims, 8 Drawing Sheets 
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REAL TIME CORRECTION OF AIRCRAFT 

FLIGHT CONFIGURATION 

ORIGIN OF THE INVENTION 

5 

This invention was made by an employee of the U.S. gov- 
ernment. The U.S. government has the right to make, use 
and/ or sell the invention described herein without payment of 
compensation therefor, including but not limited to payment 
of royalties. 10 

FIELD OF THE INVENTION 

This invention relates to real time monitoring, analysis and 
display of aircraft flight configuration variables 15 

BACKGROUND OF THE INVENTION 

An aircraft that is ascending following takeoff or descend- 
ing on approach will have measurable kinetic energy and 20 
potential energy components, and these components will 
change with time in measurable manners. Desirable energy 
states for both takeoff and landing can be determined from 
aircraft manufacturer guidance for these phases of flight. For 
example, where the approach occurs at an airport with an 25 
operable and reliable instrument landing system (ILS), the 
ILS may provide data recorded on the aircraft, to serve as a 
standard for comparing observed kinetic and potential energy 
components for an aircraft near the ground, below 2500 feet 
altitude and for an assumed straight glide slope path to a 30 
touchdown site. If the airport has no operable and reliable 
ILS, or if the aircraft is not near the ground, another mecha- 
nism for providing a standard for measurements or estimates 
is needed. On takeoff, where no electronic guidance compa- 
rable to a guide slope sensor is available, the aircraft climb 35 
profile can be compared to manufacturer guidance or to 
observed performance for recorded aircraft departures from 
the particular airport. 

The airline industry has become concerned with the prob- 
lem of unstable aircraft approaches, because approach and 40 
landing accidents often begin as unstable approaches. An 
“unstable approach” is often defined as an approach where, 
below a threshold altitude (1000 feet for IFR and 500 feet for 
VFR), the aircraft is not established on a proper glide slope 
path and with a proper air speed, with a stable descent rate and 45 
engine power setting, and with a proper landing configuration 
(landing gear and flaps extended). Airlines have developed 
approach procedures that call for abandonment of an 
approach that is determined to be unstable. 

Development and testing of methods for detecting atypical 50 
flights by N.A.S.A. has revealed that high aircraft energy 
during an arrival phase (below 10,000 feet but before begin- 
ning an approach) is the most commonreason for a flight to be 
identified as atypical or out of a statistically normal range. An 
atypical high energy arrival phase often corresponds to air- 55 
craft kinetic energy and/or potential energy that requires dis- 
sipation of 10-30 percent more energy than is required for a 
normal flight arrival phase. A normal arrival phase may cor- 
respond to about a 3-1 slope guide path and decelerating to an 
airspeed of about 250 knots during descent through 10,000 60 
feet altitude to a standard reference speed around 2,500 feet 
altitude, when an approach begins. 

More than half of the high energy arrivals identified by 
atypicality analysis were brought under control within stabi- 
lized approach criteria; some of the remainder of the high 65 
energy arrivals were abandoned. In contrast, where these 
findings were used to define and search for a high-energy 
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arrival exceedance, about three times as many excedances 
were detected (after touchdown); and the resulting unstable 
approaches were found to occur more frequently than the 
recoveries. 

It may be possible to identify a first class of high energy 
arrivals where recovery and subsequent stabilization is pos- 
sible and relatively easy, and to identify a second class of high 
energy arrivals in which recovery and subsequent stabiliza- 
tion is likely to be difficult or impossible. However, the 
present procedures for determining presence of a normal 
approach include an electronic glide slope that extends lin- 
early from the end of a target runway to the aircraft, whereas 
a typical (normal) aircraft approach path is curved and fol- 
lows the electronic glide slope only from about 1,800 feet 
above the field to the end of the runway. 

A 3-to-l glide path slope, corresponding to decrease of 
1,000 feet in altitude for every 3 nautical miles horizontal 
travel, is often desirable during an arrival phase. The air speed 
is 250 knots or less by regulation below 10,000 feet, and the 
aircraft decelerates to a reference speed before joining the 
approach path. These parameters are directly available but are 
unlikely to prove to be the only relevant parameters in deter- 
mining whether a flight arrival phase is normal or other than 
normal. 

Presently, analysis of the possibility of a normal aircraft 
approach or other normal aircraft maneuver is provided after 
the fact, in an off-line comparison of the present aircraft flight 
parameter value(s) with a collection of corresponding refer- 
ence flight parameter values that provide a “normal” aircraft 
maneuver. This off-line approach might be improved upon by 
reference to charts and tables, provided in the aircraft cabin; 
but the time required for such comparisons is realistically 
measured in minutes, and this does not allow for a real time 
response. 

What is needed is a system for monitoring, analyzing and 
displaying, in real time, a present flight configuration, of an 
aircraft and comparison of this configuration with one or 
more reference flight configurations for the same or similar 
circumstances. Preferably, the system should include the 
effects of aircraft performance constraints and should indi- 
cate whether a reference flight configuration can be achieved 
within a selected time interval, which may vary from one 
flight parameter to another and from one flight phase to 
another. 

SUMMARY OF THE INVENTION 

These needs are met by the invention, which provides 
measurements or estimates of present flight parameter (FP) 
values and compares the present value(s) with corresponding 
reference flight parameter values for similar circumstances. 
For a given flight parameter and a reference (ideal) trace or 
reference band for that parameter, a recovery band is deter- 
mined and displayed, in real time, indicating whether the 
present flight parameter value can recover to a reference trace 
or band value within a selected recovery time interval (nor- 
mally 2-15 sec), including the effects of aircraft system per- 
formance constraints on the rate of recovery. The aircraft 
performance constraints are optionally expressed in terms of 
one or more constrained differential equations involving rel- 
evant time derivatives of flight parameter values. Where 
recovery to a reference flight parameter value is not possible 
within the recovery time interval, the system optionally pro- 
vides one or more alternative responses for the pilot in com- 
mand in order to bring subsequent flight parameter values 
within the recovery band RB. The recovery band can be 
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generated and provided before the RB is used or can be 
generated at the time of its use. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 

FIGS. 1A and IB illustrates an environment in which the 
invention can be practiced. 

FIG. 2 graphically illustrates drag force developed versus 
percentage activation of an aircraft drag appliance. 

FIGS. 3 and 4 graphically illustrate a screen displaying to 
present FP values, a trace and/or band of reference FP values, 
and corresponding recovery band traces determined accord- 
ing to the invention. 

FIGS. 5A/5B graphically illustrate recovery band bound- 
ary trace determination for a reference trace and a reference 15 
band. 

FIGS. 6 and 7 are flow charts of procedures for practicing 
the invention. 

FIGS. 8-16 graphically illustrate examples of constrained 
FP values u and the applicable constraints for various con- 20 
straint classes. 

DESCRIPTION OF BEST MODES OF THE 
INVENTION 

25 

FIGS. 1A and IB illustrate environments for an ascending 
aircraft (1A) and for a descending aircraft (IB) where the 
invention can be practiced. In FIG. 1A, an aircraft 11 is 
ascending, either after takeoff or in moving from a first flight 
altitude to a second flight altitude. In FIG. IB, an aircraft 11 30 
is descending, either in preparation for a subsequent touch- 
down or in moving from a first flight altitude to a second flight 
altitude. The aircraft 11 has at least one of an associated 
kinetic energy KE(t), potential energy PE(t), energy compo- 
nent E(t)=dl:KE(t)+d2:PE(t), where (dl,d2) are selected 35 
non-negative real numbers; energy component time deriva- 
tive (d/dt)E(t), thrust power, vertical speed, ground air speed, 
aircraft mass, height above ground, flap position, speed brake 
position, landing gear position, other drag indices, roll, pitch 
and/or yaw angles, angle of attack, or another relevant flight 40 
parameter. 

The aircraft kinetic energy and potential energy compo- 
nents are, respectively, 


KE(t)=m(tyv(t) 2 /2+oyI-oi(f)/2, 

( 1 ) 

PE(i)=m{t)-g-h(t), 

( 2 ) 


where m(t) is the instantaneous aircraft mass (taking account 
of fuel consumption), I(t) is an instantaneous moment of 
inertia tensor for the aircraft, oo(t) is an aircraft rotation vector, 
computed with reference to a center of gravity or other 
selected location determined on the aircraft, v(t)=dx/dt is the 
instantaneous aircraft velocity, g is local gravity coefficient, 
and h(t) is the instantaneous height of aircraft eg above a local 
reference height, such as local ground height. The rotational 
component of kinetic energy may be negligible or may be 
ignored for other reasons. The aircraft potential energy may 
be taken as in Eq. (2). 

An aircraft system measures or estimates or otherwise 6Q 
provides a value (referred to as a “measured value” for con- 
venience herein) 

E(t n )=d\ -KE{t n )+d2-PE(t n ) ( 3 ) 

of an energy component of an aircraft during an ascent phase, 65 
a cruise phase or a descent phase of a flight, at each of a first 
sequence of times (n=l , . . . , Nl; N1^2), where dl and d2 are 
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selected real values, not both 0; for example, the pair (dl,d2) 
=(1,0), (0,1), (1,1) or (1,-1). Optionally, the system also 
provides an energy component time derivative 

(d/dt)E(t n )=d3 ■ (d/dt)KE(t n )+d4- (d/dt)PE(t n ) ( 4 ) 

where d3 and d4 are selected real values, not both 0. The 
system also provides a reference value E(f „;ref) of the energy 
component and a reference value time derivative (d/dt) E(t' w ; 
ref) for a time, t=t'„, determined with reference to the time t w 
(n=l, . . . , Nl). The time sequence {f„} may substantially 
coincide with the sequence {t w }, or each time value t' w may be 
displaced by a calculable or measurable amount from the 
corresponding time value t w . 

Each aircraft has an associated group of drag indices, one 
for each activatable drag appliance (landing gear, wing flap, 
aileron, vertical stabilizer, horizontal stabilizer, etc.). Each 
drag index has a maximum value, where the drag appliance is 
fully activated, and has a spectrum of drag values extending 
from zero activation through partial activation to full activa- 
tion of the appliance, as illustrated schematically in FIG. 2. 
With the drag appliance completely inactivated, the corre- 
sponding drag index is normally 0. The drag force associated 
with one drag appliance is assumed to be independent of the 
drag force associated with another drag appliance, in a first 
approximation. In an approach to touchdown, for example, 
where a relatively small mount of additional drag force may 
be required for adjustment, one or more drag appliances can 
be partly or fully activated to provide this additional drag 
force, relying on information illustrated in FIG. 2 for each 
drag appliance. If the amount of additional drag force needed 
for the adjustment is greater than the maximum cumulative 
drag force associated with all the activated drag appliances, 
the aircraft will need to use additional procedures to provide 
the additional drag force, or the present flight configuration 
should be (or should have been) terminated and reconfigured. 
In practice, a drag appliance is normally inactivated or fully 
activated, for some drag appliances (e.g., landing gear), and 
has a near-continuous range of settings for other drag appli- 
ances (e.g., speed brake). The sum of the drag indices for all 
(activated) drag appliances is determined and provided as a 
supplement to the drag force(s) provided by the other aircraft 
components (fuselage, wings, etc.) 

Monitoring of thrust power developed by the engine(s) of 
the aircraft is straightforward and is an important control 
variable in change of the energy components (E(t) and (d/dt) 
E(t) defined in Eqs. (3) and (4). Thrust developed can be 
estimated using measured fuel flow rate and other relevant 
variables. 

Roll angle, pitch angle and/or yaw angle of the aircraft can 
be measured or otherwise provided using one or more gyro- 
scope mechanisms that sense present orientation relative to a 
reference orientation. Angle of attack of the aircraft can be 
measured, made available and recorded during aircraft 
motion. 

FIG. 3 is an example of a screen that facilitates practice of 
an embodiment of the invention for a selected flight param- 
eter FP(t„), such as the energy component E(t„)=dl*KE(t„+ 
d2*PE(t„), the time rate of change of the eneigy component 
(d/dt)E(t„), applied thrust, a selected drag index D(t„), and 
roll, pitch, yaw and/or angle of attack for the aircraft body. 
The screen in FIG. 3 displays first and second curves, xl and 
x2, or recovery band boundary traces FP(t w ;RB;l or 2) repre- 
senting a boundary trace or traces of a recovery band (RB), as 
a function of time, where each trace can be a continuous curve 
or can be discrete. FIG. 3 also displays a third trace x3 (or 
group of traces) representing a sequence of reference values 
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FP(t„;ref) within the recovery band RB; and a fourth trace x4 
representing the measured present FP values FP(t„), which 
may lie inside or outside the recovery band RB. The reference 
FP values trace x3 is a sequence of FP values that correspond 
to a desirable aircraft flight configuration for an environment 
similar to the present flight. Alternatively, the trace t 3 can be 
replaced by a central band of reference FP values, denoted 
FP(t w ;ref;CB), corresponding to a collection of desirable FP 
values for a given time, t=t w , as illustrated in FIG. 4. 

The trace x3 in FIG. 3 refers to reference flight parameter 
values FP(t„;ref) for a reference flight, against which the 
corresponding present flight parameter values (trace x4) is to 
be compared. It is preferable to present the RB traces, xl and 
x2, in the physical units for that flight parameter. However, all 
traces can be presented in terms of dimensionless quantities 
with a common divisor. 

The RB traces, xl and x2, are defined with reference to the 
reference FP values FP(t w ;ref), as discussed in the following. 
The RB traces, xl andx2, may be, but need not be, symmetri- 
cally arranged with respect to the reference FP values. 
Optionally, the reference FP value trace is replaced by the 
central band FP(t w ;ref;CB) and has central band boundary 
traces x3; and the RB boundary traces, xl and x2, can be 
referenced to the respective central band boundary traces, 
x3A and x3B, as illustrated in FIG. 4. 

For purposes of illustration, each trace shown in FIGS. 3 
and 4 is assumed to be discrete, for example, corresponding to 
the sequence of measurement times t=t„, Optionally, the 
present FP values FP(t w ) are provided (measured or esti- 
mated) in real time as the aircraft proceeds. The values FP(t„; 
ref), FP(t w ;ref;CB;3A and 3B) and FP(t w ;RB) are preferably 
determined beforehand and are available in an on-board or 
off-board computer. 

The RB boundary traces for a particular FP are defined with 
reference to a recovery time interval length At(FP;n;rec) as 
illustrated in FIGS. 3 and 4 and optionally depending upon 
the flight parameter, the measurement time t w , and other vari- 
ables .Fora given measurement time t w and a present FP value, 
FP(t w ), that is spaced apart from the reference FP value FP(t w ; 
ref), the present FP value FP(t w ) can change at a maximum 
time rate of change, {dFP(t=t w )/dt} max , depending upon air- 
craft system constraints, so that FP(t„.) (n’>n and increasing) 
approaches and crosses the reference FP values trace x3 at a 
time t zc (n) that will depend upon the values FP(t„) and FP(t„; 
ref), the nature of the FP itself, and any constraints imposed 
upon the FP and upon the FP rate of change. 

Appendix 1 sets forth a determination of a present FP 
value, denoted u(t;t„)=FP(t), as a function of time, beginning 
at u(t=t w ; t w ), for all constraints that can be expressed as a 
linear sum of zeroth, first and second time derivatives of the 
variable u. 

B '^c2{d 2 u/dt 1 }- , t-cl {du/dt}+c®-u1=kB" (5) 

where cO, cl and c2 are selected real numbers, at least one 
being non-zero, and B’ and B" are constraint values (at least 
one being non-zero) associated with that FP value and with 
that class of aircraft. An FP constraint can be represented by 
any reasonable expression, including but not limited to the 
expressions set forth in Appendix 1 . 

As a first example, Eq. (5) with cl=c2=0 and c0>0, corre- 
sponds to a situation in which the variable u is limited by an 
upper bound B" and/or by a lower bound B’, preferably with 
B , <0<B M , with no explicit limit on rate of change; an example 
of this constrained parameter value is aircraft flap angle (for 
drag control) or aileron angle (for control of turn radius). A 
second example is Eq. (5) with c0=c2=0 and cl>0, which 
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corresponds to a limitation on time rate of change of a vari- 
able such as angle of attack. A third example is Eq. (5) with 
cl=0 and c0>0 and c2>0, which may correspond to a limita- 
tion on radial acceleration of u, such as turn rate for the 
5 aircraft. Appendix 1 sets forth 13 classes of “constrained 
solutions” u(t;t„) of the constraint relations (5) for t^t w , 
where the initial value u(t=t w ;t w )=FP(t w ), the initial rate of 
change {du/dt}^ and/or the initial acceleration {d 2 u/dt 2 } m is 
provided, where needed, in connection with Eq. (5). Altema- 
1 0 tively, one or more constraints on aircraft performance can be 
set forth as individual limits on selected flight parameters 
Where u(t=t„;t„),=FP(t w ) differs from FP(t„;ref), a con- 
strained solution u(t;t„) is chosen that causes the solution 
u(t;t w ) to approach FP(t w ;ref) (preferably monotonically) for 
15 t>t„: 

{u(t;t„)-FP(t„; ref)} (d/dt){u(t;t n )-FP(t„; ref)}<0. (6) 

The quantity 

2Q &u(t n )=u(t n ;t n )-FP(t n ;ref), (7) 

with if >n and increasing, is analyzed, and a first zero-crossing 
value, t=t zc (n)>t„ is estimated for which Au(t=t zc (n))=0. 
Determination of the zero crossing time t zc (n) is straightfor- 
ward, because FP(t„;ref) is known or specified and u(t;t w ) (for 
25 ti^t w ) is determinable from the results in Appendix 1 or other 
relevant constraint expressions. The system then determines 
if a specified recovery time interval length At(FP;rec) (depen- 
dent upon the present phase of the flight) satisfies 

30 Atl(t n )={t zc (ri)-t n }-At(FP;rec)^0. (8) 

The recovery time interval length may vary with the particular 
choice of FP and with the phase of the flight, but often lies in 
a range 2-15 sec. If the inequality in Eq. (8) is satisfied, the 
35 initial value u(t;t„)=FP(t w ) lies within the recovery band RB, 
and the present FP value can recover to a subsequent refer- 
ence FP value FP(t„,,;ref) (n’>n) within the recovery time 
interval. If the inequality in Eq. (8) is not satisfied, the initial 
value u(t=t w ;t„)=FP(t„) does not lie within the recovery band 
4Q RB, and additional time is likely required for recovery. The 
algebraically largest (or smallest) constraint equation solu- 
tion initial value, u(t;upper) or u(t;lower) for t=t w , for which 
the inequality in Eq. (8) is satisfied determines one value 
uO(t=t w ;t„) on the recovery band boundary RB (tracexl orx2). 
45 The locus of all points on one or the other of the boundary 
traces xl or x2 defines that recovery band boundary (RBB) 
trace. 

Determination of Recovery Band Boundaries 

50 As noted in Appendix 1 , to determine the coordinates of a 
point on the recovery band boundary (RBB), it is more con- 
venient and straightforward to begin at the coordinates (t,u(t)) 
of a point on the FP reference curve FP(t;ref) and to proceed 
55 backward in time, using a modified time variable x, such as 

T=t zc {n)+t n -t, (9) 

to determine a boundary point on the RBB. Beginning at a 
point on a flight parameter graph with initial coordinates 

60 {t,u{t))={t z Jn\ FP{t zc {n); ref)) (10) 

on the FP reference curve, a constrained solution is run in 
reverse to a point having final coordinates 

(t,u(t))=(t zc (n)-At(FP;rec), u(t zc (ri)At(FP;rec))), (11) 

65 

which is identified as a point on the recovery band boundary 
RBB. This process is repeated for a sequence of initial coor- 
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dinates (t zc (n ), FP(t zc (n);ref) on the reference curve, prefer- 
ably (but not necessarily) for an “initial” temporal value 

t=U»), ( 12 ) 

corresponding to an initial time t=t„ for a constrained solution 
u(t;t w ). That is to say, the choice for the value t zc (n) may be 
arbitrary but preferably corresponds to one of the sequence of 
the measurement values, t=t w , for the FP. 

A recovery band boundary trace is thus expressible as a 
sequence of coordinate pairs {(t zc (n)-At(FP;rec), u(t zc (n)-At 
(FP;rec))}, where the time value t zc (n) is arbitrary. The choice 
of the temporal value t zc (n) determines coordinates (t=t zc (n)- 
At(FP;rec), u(t=t zc (n)-At(FP;rec)'t w ) of the RBB point corre- 
sponding to that point on the FP reference curve. As noted 
above, the upper and lower RBB traces, xl and t 2, are not 
necessarily symmetrically arranged relative to the reference 
trace FP(t„;ref))} so that these two traces may need to be 
determined separately. FIG 5A illustrates this procedure 
where a single FP reference trace FP(t;ref) is used. 

FIG. 5B illustrates use of a reference band FP(t„;ref;CB) 
with upper and lower boundaries, rather than a reference trace 
FP(t w ;ref;CB) to determine the recovery band boundaries. 
The upper trace RBB(U) or lower trace RBB(L) of the RBB 
is determined only with reference to the upper RBB or lower 
RBB, respectively, as indicated in FIG. 5B. 

FIG. 6 is a flow chart of a procedure for determining an 
upper (or lower) RBB, which can be determined before a 
flight or during a flight. In step 61, a sequence of two or more 
FP reference values FP(t w ;ref) is provided, with correspond- 
ing temporal values, t=t„ (n=nl, n2), and a recovery time 
interval length At(FP;rec) is provided, for a selected flight 
parameter FP. In step 62, a constrained solution u(t;t„.) for the 
FP is determined or otherwise provided for a sequence of two 
or more initial temporal values, t=t w ., having at least one initial 
value, u(t=t w ,;t w ) that is displaced from a corresponding ref- 
erence value FP(t w .;ref), either through the considerations in 
Appendix 1 or by another suitable prescription of constraint 
(s) on the FP. In step 63, a corresponding reference curve- 
crossing time value, t=t zc (n), is determined where the con- 
strained solution u(t;t w ) crosses a reference curve FP(t;ref), or 
crosses a selected boundary trace of a reference band FP(t; 
ref;CB), for each of at least two of the sequence of constrained 
solutions u(t;t w ). In step 64, beginning with a point on the 
reference curve or reference band boundary trace with coor- 
dinates (t zc (n) u(t zc (n);t w )), the system computes the coordi- 
nate pair (t zc (n)-At(FP;rec), u(t zc (n)-At(FP;rec);tJ) for at 
least two time values, t=t w (n=nl, n2), and identifies this 
coordinate pair with a point on the recovery band boundary. 
This may be repeated for each constrained solution {u(t;t„)}„ 
in the sequence to determine a segment of the RBB in step 65 
(optional). 

Where the reference FP values FP(t m ;ref) are replaced by a 
central band of values FP(t w ;ref;CB), the quantity FP(t m ;ref) 
in steps 61, 62 and 63 is replaced by an upper reference band 
boundary trace FP(t;ref;CB;3A) or by a lower reference band 
boundary trace FP(t;ref;CB;3B) for the central band in the 
flow chart in FIG. 6, and the invention can be applied to a 
central band FP(t;ref;CB) as well as to a single reference trace 
FP(t;ref). 

This procedure is preferably implemented beforehand and 
provided as part of a software package in a computer (on- 
board or offboard) that is used on a flight. However, this 
procedure may also be implemented during the flight, in real 
time or near real time. This procedure would be implemented 
separately for each flight parameter FP of interest, and each 
FP will have its own characteristic constraint solutions, which 
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may vary with the phase of the flight, with the initial value 
u(t=t w ;t w ), and with the local environment (origin airport, 
waypoint, destination airport, etc.). 

5 Use of the Recovery Band 

When the recovery band and its boundaries (for a given 
phase and a given flight parameter FP) have been determined 
or otherwise provided, an aircraft in flight or proceeding 
io toward take off may use the RB, as indicated in a flow chart in 
FIG. 7, to determine if recovery to a reference curve or ref- 
erence band is possible. In step 71, a recovery time interval 
length At(FP;rec) and corresponding recovery band boundary 
traces FP(t;RB;l) and/or FP(t;RB;2) are determined or oth- 
1 5 erwise provided, for a selected flight parameter FP and for an 
appropriate present phase of a flight. In step 72, a present 
flight parameter value FP(t m ) of the flight parameter FP is 
provided. In step 73, the system determines if FP(t m ) is within 
or on a boundary of the recovery band, for example, 

20 FP(t„;RB; 1 ) ^FP(t m )^FP(t„;RB;2)). ( 13 ) 

If the answer to the query in step 73 is “yes,” the system 
interprets this condition as indicating that the flight parameter 
FP can (or is likely to be able to) recover to the reference curve 
25 or reference band within a recovery time interval length, 
beginning at the present time, t«t m , in step 74. If the answer to 
the query in step 73 is “no,” the system interprets this condi- 
tion as indicating that the flight parameter FP is unlikely to be 
able to recover to the reference curve or reference band within 
30 a recovery time interval length, beginning at the present time, 
t»t m , in step 75. In step 76, the system optionally recommends 
at least one action to subsequently bring the flight parameter 
value within the recovery band. The system then proceeds to 
step 77 (optional) and displays a selected section of the recov- 
35 ery band and one or more contemporaneous value(s) FP(t m ), 
with reference to the recovery band. The procedures in the 
flow charts in FIGS. 6 and 7 can be performed together. 

This procedure is implemented din real time, or near-real 
time, as the flight progresses. Where recovery of K flight 
40 parameter values FP^(t m ) (k=l, . . . , K;Ki^2) is a concern, 
each flight parameter will have its own recovery time interval 
length At(FP^;rec) and its own recovery band boundary 
traces, FP^(t;RB;l) andFP£(t;RB;2); andjoint recovery of the 
K flight parameters should be examined with reference to the 
45 K recovery bands. 

In order to generate a continuous trace of a reference curve 
or reference band boundary, such as FP(t;ref) or FP(t;ref; 
CB;1 or 2), or a continuous trace of a recovery band boundary, 
such as FP(t;RB;l or 2), one may need to extend the under- 
50 lying definition of the reference or recovery function from 
discrete points to a continuum (here, of temporal values t). 
Given N function values f(t„) (n=l, 2, . . . , N; N^2), with 
t x <t 2 < . . . <t N for definiteness, the estimation function 

NN 

F(t;est)=llfitJ%{(t-t n )/(t m -t n ) } ( 14 ) 

m=l n=l, n^m 

is continuous in the range t 1 and beyond and satisfies 

60 F(t=t k ;est)=f(t k ). ( 15 ) 

Because the measurement times are spaced apart, none of the 
denominators in F(t;est) is 0, and the magnitude of F(t;est) is 
bounded. The definition for F(t;est) can be extended to any 
65 finite number N of time points. 

More generally, one can define an estimation function F * (t; 
est) as a sum of two or more continuous characteristic func- 
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tions F*(t;k) (k=l, . . . , N; N^2), where the characteristic 
function component F*(t;k) satisfies 

F*{t=t p ;k)=At k ) (p=k)= 0 (p+k). (16) 

APPENDIX 1 

Some Constraint Equations for Aircraft Performance 


At=t-t n (^0), (A-5c) 

with uO unconstrained. FIG. 9 illustrates the upper limit and 
lower limit solutions, u(t;upper) and u(t;lower), and a permis- 
sible constrained solution u(t), corresponding to Eqs. (A-5). 
Optionally, one of the constraints can be deleted. The con- 
strained solution u(t;t„) for t^t w can be implemented and 
analytically using the functional 


Constraints on aircraft performance can arise from limita- 
tions in engine, airfoil and system performance, limitations 
imposed by an airline on its own flights, limitations imposed 
by a regulatory agency (such as the F.A.A. or the N.T.S.B.), 
and other limitations. For purposes of this analysis, a con- 
straint on a flight parameter FP is assumed to be expressible in 
terms of zeroth, first and second order linear differential equa- 
tions in a constrained solution u(t;t„) for the flight parameter 


FP: 

B^c2(d 2 u/dP)+cl (du/dt)+cO-u^B", (A-l) 

u{t=t n ;t n )=uO(t n )=FP(t n ) , (A-2a) 

{d/dt)u(t=t n ;t^)=u 'O(^) (optional; dependent on t„), (A-2b) 

{cP/dt 2 )u(t-t n ;t n )=u"0(t n ) (optional; dependent on t„), (A-2c) 


where c2, cl and cO are selected real numbers, at least one of 
which is non-zero, BO and B1 are selected real numbers. 12 
classes of constraint equation solutions are analyzed here, 
depending upon which of the coefficients c2, cl and cO are 
non- zero. For most relevant constraints (A-l), BO<0<B1. 

Class 1: cl=c2=0; cO^O. This class includes limitations on 
an aircraft flap angle (or drag control) or an aileron angle (for 
aircraft turn radius). Without loss of generality, the coefficient 
cO is assumed to have the value+1. The quantity u is con- 
strained by 

B0^u(t)^Bl, (A- 3 a) 

BO^u(t=t n ;t n )=uO^Bl, (A-3b) 

is equal to u between BO and Bl, is limited at the lower end by 
the value BO and is limited at the upper end by the value Bl 
(i^BO), as illustrated in FIG. 8. The constraints in Eq. (A-3) 
are implemented analytically by the functional and constraint 
solution 

u-^*a{u;a l;b 1 }=al +{ I u(t;t n )+bl l-lM(r;r„)-MI}/2 (A-4a) 


u->a{u;a2;b2} 

(A- 6a) 

a2(t)=(B0+Bl)(At)/2+u0, 

(A- 6 b) 

b2(t)=(Bl-B0)(At)/2+u0, 

(A- 6c) 

u(t;[ower)=B0(At)+u0^u(t;t n )=u'0- (At)+u0 ^Bl(At)+i 
u0=w(?;upper), 

(A-6d) 

The constraints set froth in Eq. (A-6d) correspond to Eq. 

(A-l). 

Class 3 : c0=cl=0,; c2^0. This class includes a limitation on 
second order rate of change of the flight parameter value, with 
uO and u’O unconstrained. Without loss of generality, the 
coefficient c2 is assumed to have the value+1 . The functional 
and constraint solution for this class are given by 

B0^u”0^Bl, 

(A- 7) 

u^*a{u;a3;b3}, 

(A- 8a) 

a3 ( t)=(B0+B 1 )■ (Af) 2 /4+w '0- (At)+u0 , 

(A- 8 b) 

b3(f)=(Bl-B0)-(At) 2 /4, 

(A- 8c) 

u(t;lower)=B0- (At) 2 /2+u '0(Af)+u0 ^ u{t;t n )=u "0 ■ (At) 2 ! 
2+u '0(Af)+«0 ^B1 • (At) 2 /2+ul -(Ar)+w0=M(r;up- 
per). 

(A-8d) 


The upper branch and lower limit solutions and a permissible 
constrained solution u(t;t w ) are illustrated in FIG. 10. 

Class 4: cO^O, cl*0, c2=0, c0/cl=a (*0). This class 
includes a combined limitation of a sum of the zeroth and first 
derivatives. The constraints are expressed as 

BO^du/dt+au^Bl, (A-9) 

and the function and constraint solution of this class are given 
by 

u^a{u;a4;b4}, (A-lOa) 


10 


15 


20 


25 


30 


35 


40 


45 


al=(B0+Bl)/2, 


(A-4b) 


a4(r)=((2?l+2?0y2a){l-exp(-aAr)}+w0-exp({-«Ar) (A- 10b) 


bl=(Bl-B0)/2, (A-4c) 

together with any other relevant initial conditions for u(t;t„) 
and its derivatives. Optionally, one of the constraints can be 
deleted, by formally setting B0=-oo or Bl=+oo. The left hand 
constraint and right hand constraint in Eq. (A-3a) would be 
applied when the quantity u is to be driven toward smaller 
values and toward larger values, respectively. 

Class 2; c0=c2=0; cl^O. This class includes a limitation on 
rate of change of angle of attack or rate of change of vertical 
or horizontal stabilizer angle. Without loss of generality, the 
coefficient cl is assumed to have the value +1 . FIG. 8, with u 
replaced by du/dt, illustrates the affect of this constraint on the 
variable du/dt. The constrained solution u satisfies. 

u(t;[ower)=B0At+u0^u(t;t n )=uX)(At)+u0 ^BlAt+ 

w0=w(r;upper), (A-5a) 


M(£)=((2?l-2?0y2a){l-exp(-aAr)} (A- 10c) 

«(r;lower)=(2?0/a){ 1 -exp(-aAr)+«0-exp{-aA?)} =u{t; 

O = ( M '0/a){ 1 -exp(-aAr)+«0-exp{-aA?)} ( Bl/a ) 
{l-exp(-aAr)+«0-exp{-aAr)=M(r;upper). (A-lOd) 

Here, uO and u’0 are constrained by Eq. (A-9). The upper 
branch and lower limit solutions and a permissible con- 
strained solution u(t;t w )are illustrated in FIG. 11. 

Class 5: c0=0, cl^O, c2*0, cl/c2=a(^0). The constraints 
are expressed as 

a 2 B0 ^kcPu/dP+a {du/dt)1=ka 2 B 1 . (A- 11a) 

The constraint relation (A-l la) can be re-expressed as 

a 2 B0 <exp(-aAr) {d/di) { exp(cit)(du/dt) } ^ a 2 B 1 , (A- lib) 


BO^uV^Bl, 


(A-5b) which has a constrained solution expressible as 
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50{aAf-(l-exp(-aAr)}+aa'0{l-exp(-aAr)}+«0 exp(- 
aAr)=w(7;lower) T u(t;t n ) TB 1 {a At- ( 1 -exp (-a At) 

}+au U{ 1 -exp(-a Ar)}=«(?;upper). (A- 1 2a). 

The functional and constraint solution for this class are given 


by 

u^>a{u;a5;b5}, (A- 12b) 

a5 ( t)=(B0+B 1 ){a At- ( 1 -exp (-a At)}/2+au U { 1 -exp (- 

a At)}, (A- 12c) 

10 

65(r)=(^l-^0){«Ar-(l-exp(-aAr))}/2. (A-12d) 

Here, uO and u'O are constrained by Eq. (A- 12a)). The upper 
branch and lower limit solutions and a permissible con- 
strained solution u(t;t w ) are illustrated in FIG. 12. 15 

Class 6: cO^O, cl=0, c2*0, c0/c2=€ 2 a 2 . (€ 2 =±1) For € 2 =-l, 
the constrains and initial conditions are written 

a 2 B0^d 2 u/d?-a 2 u^a 2 B\. (A-13a) 

u(t=t n ;t n )=u0, (A- 13b) 20 

(d/dt)u(t=t n ;t n )=u'0. (A- 13 c) 

The left hand inequality in Eq. (A- 11a) (BO only) is re- 
expressed as 25 

a 2 B0^Qxp(aAt) (d/dt){exp(-2aAt) (d/dt)(exp(aAt)u)} (A-14a) 

At=t-t n (^0). (A- 14b) 


Equation (A- 14a) can be temporally integrated once to pro- 30 
vide 

a 50{l-exp(-aAr)}exp(2aAf)+H , 0+a wO}exp(2aA?)^ 

(d/dt){u exp(aAf)}, (A-15a) 

which can be temporally integrated a second time to provide 35 

(50/2){exp(aAf)-exp(-aA*)}+(250+«0+«'0/a){exp 
(aAt)-exp(-aAt)l2}+uO exp(-aA/)=w(?;lower)^w 
(t;t„)=(B 1 /2) (exp(a A*)-exp (~aAt)+(2B 1 +w0+ 
u'0/a){exp(aAt)-exp(-aAt)}/2+u0 exp(-aAt)=u 
(fjupper), (A-15b) 40 

where the right hand inequality in Eq. (A- 13a) is included by 
symmetry arguments. The constrained solution u(t;t w ) in Eq. 

(A- 15b) can be implemented using the functional 

u^a{u;a6;b6}, (A-16a) 45 

a 6 (t)=((B0+B 1 )/4){exp («Af)-exp(-a Az)+((50+5 1 )/2) 
{l-exp(-ar)}+(«'0+a wO){exp(aAf)-exp(-aAf)}+ 
uO exp (-a At) (A- 16b) 

b 6 (r)=((5 1 -50)/4){exp («Ar)-exp(-a Az)+((51 -B0)/2) 

(l-exp(-ar)} (A-16c) 

Here, uO and u'O are constrained in accord with Eq. (A- 15b). 
The upper branch and lower limit solutions and a permissible 
constrained solution u(t) are illustrated in FIG. 13. 55 

Class 7: c0*0, cl=0, c2*0, c0/c2=€ 2 a 2 .(€ 2 =±l) For € 2 =+l, 
the constraints and initial conditions are written 

a 2 BO ^d 2 u/dt 2 +e 2 a 2 u^a 2 B 1 . (A-17a) 

u{t=t n ;t n )=u0, (A- 17b) 60 

(d/dt)u(t=t n )=u'0, (A-17c) 

The lower and upper solutions, the constrained solution, the 
functional a{u;a7;b7) and the functional parameters, a7 and 65 
b7, are determined by replacing the parameter a by +ia or by 
-ia (i 2 =-l) in Eqs. (A- 15b), (A-16a), (A-16b) and (A-16c). 


12 

{{B0l2)+{u'0+a wO)}{exp(*aAr)-exp(-*aAf)}+aO exp 
(-iaAt)=u(t\\o\ved)1=k u(t;t n ) =i{ (5 l/2)+(w'D+a 
«0)}{exp (iaAt)-exp(-iaAt)}+uO exp ( -iaAt=u(t ; 


upper), (A- 18a) 

u^a{u;al;bl}, (A- 18b) 

al{t)={{B0+B\)IA){sxp{iaAt)-Qxp{-iaAt)}+{u'0+a 

wO){exp(aAr)-exp(-aAr)}+wO exp (- aAt ) (A- 18c) 

bl(t)=((Bl-B0)l2) { exp(iaAt) - exp (-iaAt)}, (A-18d) 


Here, uO and u'O are constrained in accord with Eq. (A- 17a). 
The upper branch and lower limit solutions and a permissible 
constrained solution u(t;t w ) are illustrated in FIG. 14. 


Class 8: cl/c2=2a^0, c0/c2=xa 2 (%*+l). The defining 
equation becomes 

BOa 2 1!kSu/dt 2 +2a{du/dt)+xa 2 u1=kB\ a 2 . (A-l 9a) 

Define 

-al=±2^{(l-x 2 )fl 2 }, (A- 19 b) 

a2=a±J{(l~x 2 )a 2 }, (A- 19c) 

which are components of the roots of the quadratic equation 

A 2 +2aA+xa 2 =0, (A-20a) 

A=a{-l±J(l-c)}. (A- 20 b) 


The left hand side of the defining equation (A- 16) is re- 
expressed as 

50 a 2 ^exp{(- ( a 1 +a2)At}(d/dt){exp ( a 1 At) (d/d t) {exp 

(a2At)a}} (A-21) 

Equation (21) can be temporally integrated once to yield 

{a 2 50/(al+a2)}{exp{(al+a2)Ad-l}=exp(«lAr) 

(d/dt){exp(a2At)u}-u'0+a2 uO (A-22a) 

Equation (A-22a) can be temporally integrated a second time 
to yield 

{a 2 50/(al+«2)}{(l-exp(-a2A?))/a2-(exp(-(«2A?)- 
exp(-(al+a2)A?))/al}+(w'0+a2 wO)(exp(-a2Ar)- 
exp(-(a 1 +a2)At)!a l+wO=w(r;lower)^«(f;£ n ) 
^{a 2 51/(al+«2)}{(l-exp (-a2At)/a2-(exp(- 
(«2At) -exp (~(al+a2)At))l a\}+(u'0+a2 «0)(exp 
(-a2Ar)-exp(-(al+a2)Ar)/al+«0=«(r;upper), (A-22b) 

where the right hand inequality in Eq. (A-22b) is included by 
symmetry arguments. The constrained solution u(t;t w ) in Eq. 
(A-22b)) can be implemented using the functional 

w-*a{w;«8;68}, (A-23a) 

a8(r)={a 2 (50+51)/2(al+a2)}{l-exp(-a2At))/a2-(exp 
(-(«2Af)-exp(-(al+a2)Ar))/al }+(w'0+a2 «0)(exp 
(-«2At)-exp(-(al+a2)Ar)/al+«0 (A-23b) 

68(f)=(a 2 (51-50)/2(al+a2)}{(l-exp(-a2Af))/a2- 

(exp(-(a2Ar)-exp(-(al+a2)Ar))/al}. (A-23c) 

Here, uO and u'O are constrained in accord with Eq. (A-l 8a). 
The upper branch and lower limit solutions and a permissible 
constrained solution u(t;t w ) are illustrated in FIG. 15. 

Class 9: cl/c2=2a^0, c0/c2=xa 2 (x=+l). The defining 
equation becomes 

50a 2 ^d 2 u/dt 2 2a(du/dt)+a 2 u=L[u]^Bl a 2 , (A-24) 

which has a degenerate solution u=exp(-aAt) (a double root 
of L[u]=0). Assume that a second, independent solution of 
L[u]=0 has the form 
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u(t;t„)=w(t;t n )ex p (-a At) . (A-25) 

The solution component w(t) is found to satisfy 

d 2 w/dt 2=0, (A-26a) 

5 

w(r;r„)=uO+wtDAr (A-26b) 

The constraints in Eq. (24) can be re-expressed as 

50{{l-exp(-tfAr)}-aArexp(-aAr)}+(w0+w'0Ar)exp(- 

aAt)=u(t;iowQr)^u(t;t n ) HkBl {{ 1 -exp(-aAf)}- 10 

flArexp(-aAr)}+(wO+w'OAdexp(-flAr)=M(r;upper) (A-27) 

The constrained solution u(t;t w ) in Eq. (A-27)) can be imple- 
mented using the functional 

u^a{u;a%;b%}, (A-28a) 15 

a9(t)=(B0+B 1 ) { { 1 -exp(- a A?) }-aA?exp (-a At) }/2+ 

(wO + w OArjexp (-aAt) ,b9(t)=(B 1 -50){ { 1 -exp 
(-aAr)}-aArexp(-aAr)}/2. (A-28c) 

20 

Here, wO and w’0 are constrained in accord with Eq. (A-27). 
The upper branch and lower limit solutions and a permissible 
constrained solution u(t;t w ) are illustrated in FIG. 16. 

Class 10: BO^u^Bl and B2^du/dt^B3. This situation is 


described by two simultaneous equations, 

50 , (A-29a) 

B2At^u(t;t„)-uO^B3At, (A-29b) 

which may be implemented by the functional 3Q 

u-*a{u;a9;b9}, (A-30a) 

a 1 0(r)={max{50, 52AM-w0)+min {SI , B3At uO}}!2, (A-30b) 

b 1 0(f)={min{Sl, B3At w0}-max{S0, B2At}}/2. (A-30c) 

Class 11: BO^u^Bl and B4^d 2 u/dt 2 ^B5. This situation 
is described by two simultaneous equations 

B0%u(t;t n )^B\, (A-31a) 

B4(At) 2 l2%u{t;t n )-u'0At-u0^B5{At) 2 l2, (A-31b) 

which may be implemented by the functional 

u^a{u;al0;b\0}, (A-32a) 

al l(?)={max{S0, S4(Ar) 2 /2+w'0Ar+w0}+min{Sl, 

B5(Atf/2+uOAt+uO}/2, (A-32b) 

b 1 0(r)={min{51, B5(At) 2 /2+uOAt+uO}-ma.x{BO, 

B4(At) 2 /2+uOAt+uO)/2. (A-32c) 

Class 12: B2^du/dt^B3 and B4^d 2 u/dt 2 ^B5. This situ- 
ation is described by two simultaneous equations 

B2At^u(t;t„)-u'0^B3At, (A-33a) 

B4(At) 2 /2^u(t;t„)-uX)At-u0 ^B5(At) 2 /2, (A-33b) 

which may be implemented by the functional 

u^a{u;al\;bl\}, (A-34a) 

a 1 2(?)={max{S2+« U, S4(Ar) 2 /2+M'0Ar+«0}+min{S3+ 

uO, B5{At) 2 }l2+uOAt+uO}l2, (A-34b) 

b 1 2{t)={mm{B3+u U, S5 {At) 2 !2+u '0AM-w0}-max{S2+ 

uO, B4(At) 2 l2+uOAt+u())}l2. (A-34c) 

Class 13: BO^u^Bl and B2^du/dt^B3 and B4^d 2 u/ 


dt 2 ^B5. This situation is described by three simultaneous 
equations 
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B0^u(t;t n )^Bl, (A-35a) 

B2At^u(t;t n )-uO^B3At, (A-35b) 

B4At^ u(t;t n )-u f 0At-u0 ^B5At, (A-3 3c) 

which may be implemented by the functional 

w-*a{w;al2;M2}, (A-3 6a) 

a\3(t)={max{B0, B2At+uO, B4(At) 2 /2+uOAt+uO}+ 

min {51, B3At+uO, B5(At) 2 /2+uOAt+uO}}/2, (A-3 6b) 

b\3(t)={mm{Bl, B3At+uO, B5{At) 2 l2+u'0At+u0}- 

max{50, B2At+uO, B4{At) 2 l2+uOAt+uO}}l2. (A-3 6c) 


The defining differential equations (e.g., Eq. (A-l)) permit 
constrained solutions for t^t w as well as for X=X n . The con- 
strained solutions and recovery band boundary values FP(t; 
RB) are more conveniently used, beginning with an FP “ini- 
tial value” on the reference trace FP(t;ref) or on a reference 
boundary trace FP(t;ref;CB;l or 2) and computing the con- 
strained solution in a backward-running direction. This may 
be accomplished in at least two manners; (1) compute con- 
strained solutions U(x;t„) of Eq. (A-l ) where the time variable 
x runs backward, x=t0-t, and where tO is selected “initial time 
value”); or (2) compute a new constrained solution U, defined 


by 

U(t;t n )=u(x;t n ), (A-3 7a) 

T=t0 -t, (A- 3 7 b) 

U(t=tO;t n )=u(xO;t„)=ul , (A-3 7c) 

(dU/dt) lmt0 =-u t 1, (A-37d) 


where u(t;t„) is the constrained solution computed in the 
preceding, tO is a selected “initial time value,” and ul and u'l 
are selected initial values. Here, the approach (2) with the 


particular choices 

x=(t zc (n)+t„)-t(t„ „(«)), (A-3 8a) 

r(initial) =?.,„(«), (A-3 8 b) 

T(initial)=f M , (A-38bc) 

u 1 =FP(t= t zc (n) ;ref) , (A-38d) 

u'\=- (d/dt)FP(t=t z Jn); ref) , (A-3 8e) 


will be used, but the approach (1), mutatis mutandis, should 
yield substantially the same results, The constrained solution 
U(x) runs backward in time, beginning at the initial coordi- 
nates (x=t„, U (x=t„)=FP(t=t zc (n) ;ref) and proceeding to a 
coordinate location 

(x=t n +At(FP;rec) , U(x=t n +At(FP;rec);t n ))=(t=t zc (n)- 

At(FP;rec), u(t=t zc (n)-At(FP;rec);t„), (A-39) 

which is identified as a boundary point for the recover band 
RB. 

What is claimed is: 

1. A method for monitoring and analyzing, in real time, 
variation with time of an aircraft flight parameter, the method 
comprising: 

measuring a value FP(t=t^) of a selected flight parameter 
FP for an aircraft, for a time, t=t^ that is equal to or 
adjacent to a present time; 

providing a sequence {FP(t„;ref)}„ of FP reference values 
for at least two time values, t=t n (n=nl, n2), and an 
available recovery time interval length At(FP;rec); 
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providing at least one constrained solution u(t;t„) for the FP 
on a graph (t, u(t)), with at least one initial constrained 
solution value u(t=t„;t„) that is displaced from the cor- 
responding value FP(t=t„;ref); 
for each constrained solution u(t;t„), providing a reference 
curve-crossing time value, t=t zc (n), at which u(t;t w ) 
crosses at least one trace of the reference values FP(t;ref) 
on the graph; 

computing a coordinate pair (t,u)=(t zc (n)-At(FP;rec),u(t zc 
(n)-At(FP;rec);t w )) on the graph for the at least two time 
values, t=t„ (n=nl, n2), and identifying this coordinate 
pair with at least one point on a recovery band boundary 
trace; 

connecting at least two of the points with the coordinates 
(t zc (n)-At(FP;rec),u(t zc (n)-At(FP;rec);t„)) to provide a 
portion of a recovery band boundary trace; 
providing at least first and second non-coincident recovery 
band boundary traces, (t,FP(t;rec;l)) and (t;FP(t;rec;2)), 
that are spaced apart from each other for at least one time 
value t on the graph (t, u(t)); 

displaying the at least first and second recovery band 
boundary traces, (t,FP(t;rec;l)) and (t;FP(t;rec;2)), and 
the value FP(t=t ) of the flight parameter at the time t=t 
on the graph (t, u(t)); and 

providing an operator of the aircraft with information 
which may be used to bring the flight parameter FP 
within the recovery band defined by the first and second 
recovery band boundary traces. 

2. The method of claim 1, further comprising: 
providing a value FP(t=t^;4>) of said selected FP parameter 

in a selected phase of a flight; 

determining if the value FP(t=t„;(|)) likes between, or lies 
on, said first and second recovery band boundary 
traces; 

when the value FP^H^) lies between said first and 
second recovery band boundary traces, or lies on at 
least one of said first and second recovery band 
boundary traces, interpreting this condition as indi- 
cating that said selected flight parameter is likely to be 
able to recover form the value FP(t=t p ;(|)) to at least 
one of said FP reference values within a time interval 
of length no greater than At(FP;rec). 

3. The method of claim 1, further comprising: 
providing a value FP^H^lOO ) of said selected FP param- 
eter in a selected phase (|) of a flight; 

determining if the value FP(t=tp;(|)) lies between, or lies on, 
said first and second recovery band boundary traces; 
when the value FP(t=t p ;(|)) does not lie between said first 
and second recovery band boundary traces, and does not 
lie on at least one of said first and second recovery band 
boundary traces, interpreting this condition as indicating 
that said selected flight parameter is unlikely to be able 
to recover from the value FP^^;^) to at least one of 
said FP reference values within a time interval of length 
no greater than At(FP;rec). 

4 . The method of claim 1 , further comprising providing 
said constrained quantity u(t;t„) for said FP value as a dimen- 
sionless ratio FP(t)/FP0, where FPO is a selected measured or 
estimated value of said FP value. 

5 . The method of claim 1 , further comprising providing 
said constrained value u(t;t„) for said flight parameter FP as a 
solution to a selected constrained differential equation 

B '^cl{d 2 u/dp- }+cl {du/dt}+cO- u^B ", 
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where the coefficients cO, cl and c2 are first, second and third 
numerical constants, at least one being non-zero, and B' and 
B" are fourth and fifth numerical constants, at least one being 
non- zero. 

5 6. A method for monitoring and analyzing, in real time, 

variation with time of an aircraft flight parameter the method 
comprising: 

providing a recovery band in a graph of values (t, u(t)) of a 
measured flight parameter (FP) for an aircraft versus 
to time t, the recovery band being defined by first and 
second non-coincident recovery band boundary traces 
that are spaced apart for at least one time value, and 
providing a recovery time interval length At(FP;rec), 
having a positive temporal value, where the flight 
1 5 parameter, having a value FP(t=tp) at a time , i s likely 

to be able to recover to a reference flight parameter value 
FP(t';ref), lying in a band of reference flight parameter 
values FP(t';ref;CB), within a time interval given by 
t /? ^t , ^t p +At(FP;rec), if the flight parameter value lies 
20 between, or lies on, the first and second recovery band 
boundary traces; and 

displaying at least one of the recovery band on the graph 
and the recovery time interval length At(FP;rec). 

7. The method of claim 6, further comprising providing 
25 said reference flight parameter band as a curve having a width 

of substantially 0. 

8. The method of claim 6, further comprising providing 
said quantity w(t) in said graph as a dimensionless ratio 
FP(t)/FP0, where FPO is a selected measured or estimated 

30 value of said FP value. 

9. The method of claim 6, further comprising providing 
said value w(t) for said flight parameter FP as a solution to a 
selected constrained differential equation 

35 B '= c2 { d 2 w/dp- (+c 1 {dw/dt}+c 0 -w^B ", 

where the coefficients cO, cl and c2 are first, second and third 
numerical constants, at least one being non-zero, and B' and 
B M are fourth and fifth numerical constants, at least one being 

non-zero. 

40 

10. A method for monitoring and analyzing, in real time, 
variation with time of an aircraft flight parameter the method 
comprising: 

providing a recovery band in a graph of values (t, u(t)) of a 
45 measured flight parameter (FP) for an aircraft versus 
time t, the recovery band being defined by first and 
second non-coincident recovery band boundary traces 
that are spaced apart for at least one time value, and 
providing a recovery time interval length At(FP;rec), 
5Q having a positive temporal value, where the flight 
parameter, having a value FP(t=t^) at a time t=t p , is 
unlikely to be able to recover to a reference flight param- 
eter value FP(t';ref), lying in a band of reference flight 
parameter values FP(t';ref;CB), within a time interval 
55 given by t^f ^t^+At(FP;rec), if the flight parameter 
value does not lie between, and does not lie on, the first 
and second recovery band boundary traces; and 

displaying at least one of the recovery band on the graph 
and the recovery time interval length At(FP;rec). 

60 11 . The method of claim 10, further comprising providing 

said reference flight parameter band as a curve having a width 
of substantially 0. 

12. The method of claim 10, further comprising providing 
said quantity w(t) in said graph as a dimensionless ratio 

65 FP(t)/FP0, where FPO is a selected measured or estimated 
value of said FP value. 

13. The method of claim 10, further comprising 
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providing said value w(t) for said flight parameter FP as a 
solution to a selected constrained differential equation 

B'^c2{d 2 w/dt 1 }+c\{dw/dt}+cO-w^B", 

where the coefficients cO, cl and c2 are first, second and third 
numerical constants, at least one being non-zero, and B' and 
B" are fourth and fifth numerical constants, at least one being 
non- zero. 

14. A system for monitoring and analyzing, in real time, 
variation with time of an aircraft flight parameter, the system 
comprising a computer that is programmed: 

to measure a value FP(tp) of a selected flight parameter FP 
for an aircraft, for a time, t=t^ that is equal to or adjacent 
to a present time; 

to provide a sequence {FP(t w ;ref)}„ of FP reference values 
for at least two time values, t=t„ (n=nl, n2), and an 
available recovery time interval length At(FP;rec); 
to provide at least one constrained solution u(t;t M ) for the 
FP on a graph (t, u(t)), with at least one initial con- 
strained solution value u(t=t„;t„) that is displaced from 
the corresponding value FP(t=t w ;ref); 
for each constrained solution u(t;t„), to provide a reference 
curve-crossing time value, t=t zc (n), at which u(t;t„) 
crosses at least one trace of the reference values FP(t;ref) 
on the graph; and 

to compute a coordinate pair (t,u)=(t zc (n)-At(FP;rec),u(t zc 
(n)-At(FP;rec);t w )) on the graph for the at least two time 
values, t=t„ (n=nl, n2), and to identify this coordinate 
pair with at least one point on a recovery band boundary 
trace; 

to connect at least two of the points with the coordinates 
(t zc7 (n)-At(FP;rec),u(t zc (n)-At(FP;rec);t K )) to provide a 
portion of a recovery band boundary trace; 
to provide at least first and second non-coincident recovery 
band boundary traces, (t,FP(t;rec;l)) and (t;FP(t;rec;2)), 
that are spaced apart form each other for at least one time 
value t on the graph (t, u(t)); and 
to display the at least first and second recovery band bound- 
ary traces, (t,FP(t;rec;l)) and (t;FP(t;rec;2)), and the 
value FP(t=tp) of the flight parameter at the time t=t p on 
the graph (t, u(t)); and 

to provide an operative of the aircraft with information 
which may be used to bring the flight parameter within 
the recovery band defined by the first and second recov- 
ery band boundary traces. 

15. The system of claim 14, wherein said computer is 
further programmed: 

to provide a value FP(t=t^;(|)) of said selected FP parameter 
in a selected phase (j) of a flight; 
to determine if the value FP(t=t^;(|)) lies between, or lies on, 
said first and second recovery band boundary traces; and 
wherein the value FP(t=t^;(|)) lies between said first and 
second recovery band boundary traces, or on at least one 
of said first and second recovery band boundary traces, 
to interpret this condition as indicating that said selected 
flight parameter is likely to be able to recover form the 
value FP(t=t p ;(|)) to at least one of said FP reference 
values within a time interval of length no greater than 
At(FP;rec). 
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16. The system of claim 14, wherein said computer is 
further reprogrammed: 

to provide a value FP(t=t^,;(|)) of said selected FP parameter 
in a selected phase (|) of a flight; 

5 to determine if the value FP(t=t^; (|)) lies between, or on, 
said first and second recovery band boundary traces; 

when the value FP(t=t p ;(|)) does not lie between said first 
and second recovery band boundary traces, and does not 
lie on at least one of said first and second recovery band 
10 boundary traces, to interpret this condition as indicating 
that said selected flight parameter is unlikely to be able 
to recover from the value FP(t=t^;(|)) to at least one of 
said FP reference values within a time interval of length 
no greater than At(FP;rec). 

15 17. A system for monitoring and analyzing, in real time, 

variation with time of an aircraft flight parameter the system 
comprising a computer that is programmed: 

to provide a recovery band as a graph of values (t, u(t)) of 
a measured flight parameter (FP) for an aircraft versus 
20 time t, the recovery band being defined by first and 
second non-coincident recovery band boundary traces 
that are spaced apart at at least one time value, and to 
provide a recovery time interval length At(FP;rec), hav- 
ing a positive temporal value, where the flight parameter, 
25 having a value FP(t=tp) at a time t=t^, is likely to be able 
to recover to a reference flight parameter value FP(t',ref), 
lying in a band of reference flight parameter values 
FP(t';ref;CB), within a time interval given by t^t'^it^+ 
At(FP;rec), if the flight parameter value lies between, or 
30 on, the first and second recovery band boundary traces; 
and 

to display at least one of the recovery band on the graph and 
the recovery time interval length At(FP;rec). 

18. The system of claim 17, wherein said compute is fur- 
35 ther programmed to provide said reference flight parameter 

band as a curve having a width of substantially 0. 

19. A system for monitoring and analyzing, in real time, 
variation with time of an aircraft flight parameter the system 
comprising a computer that is programmed: 

40 to provide a recovery band as a graph of values (t, u(t)) of 
a measured flight parameter (FP) for an aircraft versus 
time t, the recovery band being defined by first and 
second non-coincident recovery band boundary traces 
that are spaced apart at at least one time value, and to 
45 provide a recovery time interval length At(FP;rec), hav- 
ing a positive temporal value, where the flight parameter, 
having a value FP(t=t^) at a time t=t^,is unlikely to be 
able to recovery to a reference flight parameter value 
FP(t';ref), lying in a band of reference flight parameter 
50 values FP(f ;ref;CB), within a time interval given by 
t^t'^t^+AftTP^ec), if the flight parameter value does 
not lie between the first and second recovery band 
boundary traces; and 

to display at least one of the recovery band on the graph and 
55 the recovery time interval length At(FP;rec). 

20. The system of claim 19, wherein said computer is 
further programmed to provide said reference flight param- 
eter band as a curve having a width of substantially 0. 



